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ABSTRACT 


Measurements  are  reported  of  the  gamma  radiation-induced  dose  in  a 
compartmented  steel  structure.  The  experimental  dose  profile  is  compared 
with  predictions  based  on  computational  techniques  in  the  Office  of  Civil 
Defense  Engineering  Guide,  as  well  as  with  computations  performed  by  the 
Monte  Carlo  technique.  Although  agreement  (within  experimental  accuracy) 
is  reached  between  the  Monte  Carlo  and  the  experimental  measurements , 
the  OCD  predictions  appear  to  overestimate  the  dose  because  of  an  incom¬ 
plete  consideration  of  the  effect  of  spacing  of  the  structure  slabs.  The  de¬ 
gree  of  conservatism  does  not  appear  to  be  serious,  however. 
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CHAPTER  1 
INTRODUCTION 

X 

The  study  of  radiation  shielding  has  generated  a  large  oody  of  literature  on  the 

effects  of  matter  in  bulk  on  a  gamma  radiation  field.  These  effects  are  normally 

♦ 

found,  both  by  theory  and  experiment,  to  be  a  reduction  of  the  dose  rate  behind 
the  matter  and  an  increase  of  the  dose  rate  in  front  of  the  matter,  i.e. ,  on  the  side 
adjacent  to  the  gamma  sources.  The  extent  of  the  effect  becomes  less  predictable 
for  the  cases  of  the  matter  forming  a  complicated  structure  such  as  one  surround¬ 
ing  the  field  point,  and  here  one  resorts  either  to  an  artful  combination  of  duct 
theories,  albedo  theories,  and  ray  theories,  or  to  the  mathematical  game  of 
Monte  Carlo  played  with  a  digital  computer. 

Such  a  problem— that  of  predicting  the  gamma  radiation  field  within  a  shield  of 
complicated  structure— has  a  practical  application.  Should  there  ever  be  a  nuclear 
attack  on  the  United  States,  radioactive  debris  from  nuclear  weapons  exploded  on 
targets  can  be  expected  to  shower  down  over  large  areas  otherwise  untouched  by 
the  effects  of  the  explosions .  In  an  idealized  theoretical  sense ,  this  fallout  can  be 
expected  to  form  a  horizontal  source  plane  over  the  ground,  perhaps  interrupted  by 
patches  suspended  by  building  roofs  and  vegetation.  These  roofs,  and  the  floors 
beneath,  may  provide  a  measure  of  protection  to  building  occupants.  Knowledge  of 
the  field  "eduction  within  the  building  is  then  of  vital  importance  in  the  evaluation 
of  fallout  shelters . 

Guides  for  the  theoretical  computation  of  such  a  field  reduction  have  been  pub- 

1  2 

lished  by  the  Office  of  Civil  Defense.  ’  These  are  based  on  empirical  extrapola¬ 
tion  of  bulk  shield  studies  of  the  radiation  transport  phenomena  and  are  designed  to 
include  such  cases  as  that  of  the  present  discussion.  The  methods  presented  in 
these  OCD  guides  are  intended  to  provide  a  close  estimate  of  the  ratio  of  the  dose 
rate  in  the  building  at  any  given  point,  to  the  dose  rate  that  would  be  experienced 


Throughout  this  report,  the  words  ”dose  rate"  refer  to  the  rate  of  creation 
of  ion  pairs  in  a  unit  volume  of  air-equivalent  medium  by  the  gamma  field .  The 
unit  of  measurement  for  dose  rate,  in  this  sense,  is  the  roentgen  per  hour. 
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in  the  absence  of  the  building  3  ft  above  the  same  field  of  fallout.  This  ratio  is 
referred  to  as  the  "reduction  factor,  "  and  its  inverse  as  the  "protection  factor.  " 
a  significant  effer*  in  this  calculation  is  that  provided  by  the  removal  of  the  por¬ 
tion  of  the  source  plane  within  the  perimeter  of  the  building  walls  to  a  position  on 
the  building  roof. 

The  purpose  of  the  present  study  was  to  measure  this  effect  experimentally  on 
an  extremely  simplified  compartmented  model  structure ,  so  that  complicating 
features  would  not  obscure  the  effect.  A  two-way  comparison  of  experiment  with 
theory  was  possible,  since,  in  addition  to  the  techniques  available  in  Refs.  1  and 
2,  a  digital  computer  calculation,  applying  a  recently  developed  program  based  on 
Monte  Carlo  techniques ,  was  made  of  the  gamma  field  for  a  central  point  in  the 
configuration  studied. 

Our  shield  configuration,  then,  was  a  structure  made  of  iron  slabs  and  re¬ 
sembling  an  egg  crate  with  three  compartments.  The  slabs  corresponded  to  the 
roof,  floors,  and  walls  of  a  simplified  building;  their  thicknesses  were  varied  to 
span  the  range  of  areal  densities  encountered  in  present-day  buildings .  A  beam 
of  radiation  from  a  strong  cobalt-60  source  was  projected  onto  a  square  cross 
section  subtended  by  the  "roof"  of  the  structure.  Data  for  various  structures, 
positions  in  the  structure,  and  angles  of  the  beam  were  obtained  by  timed  ex¬ 
posures  of  200  mr  and  10  mr  ionization  chambers;  some  of  the  more  critical  meas¬ 
urements  were  made  with  a  highly  sensitive  rate  chamber  in  a  Wheatstone  bridge 
circuit.  A  more  complete  description  of  the  experiment  is  given  in  Chapter  2. 

The  data  obtained  in  these  experiments  are  collected  in  Chapter  3 ,  with  an  ex¬ 
planation  of  the  methods  used  in  correcting  for  various  exnerimental  complications . 

In  Chapter  4,  the  experimental  data  are  analyzed  and  their  computational  tech- 

2 

niques  compared  with  those  used  in  the  OCD  Engineering  Guide .  The  analysis  con¬ 
tinues  in  Chapter  5  with  a  comparison  of  experimental  and  Monte  Carlo-calculated 
results .  Here  we  see  enough  agreement  to  justify  our  confidence  in  the  results 
derived  by  the  Monte  Carlo  technique.  Finally,  conclusions  and  recommendations 
gained  from  this  study  are  summarized  in  Chapter  6 . 
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CHAPTER  2 


DESCRIPTION  OF  EXPERIMENT 
SITE 

In  the  sprinji  of  1960,  Technical  Operations  Research  constructed  a  facility  to 
conduct  an  experimental  radiation  shielding  program  for  the  Office  of  Civil  Defense. 
The  following  year,  a  vinyl-coated  nylon  hemisphere  supported  by  a  differential  air 
pressure  was  erected  on  the  facility  test  pad  to  provide  shelter  for  all-weather 
operations.  (The  p'esent  experiments  were  performed  during  winter  months.)  As 
the  facility  has  been  discussed  in  an  earlier  report,  it  will  not  be  further  dis¬ 
cussed  here  except  for  some  details  that  are  pertinent  to  the  present  series  of 
experiments . 

Normally,  the  air  pressure  in  the  hemisphere,  which  is  maintained  by  air 
blowers,  is  14.75  psi.  The  mean  free  path  of  cobalt-60  radiation  (used  in  the  pres¬ 
ent  experiment)  is  445  ft  in  air  at  standard  atmospheric  pressure.  The  increased 
pressure  in  our  facility  slightly  reduces  this  to  443  ft.  The  dimensions  of  the  model 
structure  used  in  the  present  experiment  are  so  much  smaller  than  these  vahies, 
however,  that  one  could  reasonably  assume  no  attenuation  by  air  without  introducing 
significant  error. 

The  test  pad  of  the  facility  has  a  flat  surface  of  asphalt.  To  minimize  the  in¬ 
terference  of  gamma  scattering  from  this  surface,  the  experiment  was  conducted 
in  a  plane  approximately  4  ft  above  the  pad.  Further,  the  use  of  thick  shielding 
under  the  source  prevented  it  from  radiating  directly  to  the  floor,  thus  making  the 
effects  of  floor  scattering  negligible  in  this  experiment . 

SOURCE 

An  encapsulated  source  of  cobalt-60  was  used  for  all  experiments  conducted. 
This  source  was  compared  on  December  7,  1959,  with  one  calibrated  at  the 
National  Bureau  of  Standards  against  standard  sources  and  was  found  at  that  time 
to  have  a  strength  of  102  curies.  Based  on  a  cobalt  halflife  of  5.28  yr,  calculation 
shows  that  on  October  17,  1963,  when  these  experiments  were  initiated,  the  source 
strength  was  reduced  to  61.7  caries. 
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The  source  was  housed  in  a  Tech/Ops  Model  528  Gamma  Ray  Projector 
(Figure  1),  which  is  essentially  a  container  with  12-in.  lead  walls  and  a  curved 
duct  through  which  the  source  travels.  The  source  is  pushed  out  of  the  duct  by 
a  remotely  controlled  cable  attached  to  the  source  and  operated  by  a  crank  at  the 
control  console.  Mechanical  switches  are  built  into  the  device  and  connected  to 
indicator  llg.its  on  the  panel,  so  that  the  position  of  the  source  may  be  ascertained 
at  a  glance . 


Figure  1.  Gamma  Ray  Projector  Housing  Cobalt-60  Source 
Used  in  Experiments 
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The  gamma  ray  projector  was  seated  on  a  platform  of  concrete  blocks.  This 
platform  also  provided  space  for  a  lead  collimator,  which  was  placed  directly  in 
front  of  the  projector  aperture.  A  1/2  in.  steel  plate  on  the  platform  made  a  firm 
surface  for  projector  and  collimator. 

COLLIMATORS 

When  the  source  was  extended  from  the  projector,  it  entered  a  lead  chamber 
having  a  small  opening  that  allowed  a  beam  of  nearly  parallel  radiation  to  be  emitted. 
A  lead  slab  with  a  square  aperture  through  its  center  was  the  first  collimator.  This 
aperture  tapered  from  an  outer  2  in.  square  hole  to  an  inner  1-5/8  in.  square  hole 
and  collimated  the  beam  so  that  it  had  a  2  ft  square  cross  section  at  the  e:i'q)erimental 
target.  Figure  2  shows  the  source  assembly  with  collimator  and  with  a  cylindrical 
disk  of  lead  that  may  be  rolled  in  front  of  the  aperture.  Figure  3  is  a  schematic  dia¬ 
gram  of  the  collimator  apparatus.  Lead  and  concrete  shielding  on  all  sides  of  the 
assembly  reduced  the  surface  dose  rate  to  2  mr/hr  with  source  extended. 


Figure  2.  Projector  Housed  in  Lead  Collimator 
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A  penumbral  collimator  was  placed  10  ft  from  the  first  collimator .  This  col¬ 
limator  was  constructed  of  4  in .  thick  lead  bricks ,  the  top  layer  of  which  was  sup¬ 
ported  by  a  steel  plate,  and  stood  on  a  concrete  block  pedestal.  Since  the  beam 
widened  with  distance  from  the  point  source,  the  opening  in  this  collimator  was 
12  in.  square.  The  major  purpose  of  this  construction  was  to  cut  off  those  portions 
of  the  beam  partially  obstructed  by  the  first  collimator,  thereby  "shaping"  a  beam 
of  almost  constant  intensity;  some  small  angle  air- scattering  was  probably  also 
thus  eliminated. 

SIMULATED  BUILDING  STRUCTURE  (EGG  CRATE) 

The  projector  and  collimator  allowed  a  gamma  beam  of  a  square  cross  section 
with  a  maximum  5®  spread  to  strike  the  experimental  structure.  During  the  course 
of  the  experiments,  this  structure  was  dubbed  the  "egg  crate.  "  The  basic  unit  of 
the  egg  crate  was  four  steel  slabs,  1/2  in.  thick  and  2  ft  square,  each  supported  in 
a  vertical  plane  by  two  steel  legs.  2  in.  x  1/2  in.  x  36  in.  The  legs  were  connected 
by  a  Cl  oss  strap  and  a  base  platform  under  the  slabs ,  and  the  slabs  positioned  1  ft 
apart.  All  these  connections  were  made  by  spot- welding. 

During  the  experiments,  a  bottom  plate  section  was  added  to  this  basic  unit. 

It  was  composed  of  three  1/2  in.  thick  plates,  each  of  which  fitted  into  the  space 
between  two  adjacent  slabs  to  form  a  floor.  Side  plates  of  the  same  thickness  were 
later  added,  and  these  rested  on  tabs  projecting  from  the  legs.  Finally,  a  top 
plate  was  laid  over  the  basic  unit  to  form  a  completely  enclosed  box  with  two  in¬ 
ternal  partitions.  All  surfaces  of  the  egg  crate  were  painted  with  a  thin  coat  of 
rust-resistant  paint.  Figure  4  shows  a  view  of  the  instrumented  egg  crate  just 
prior  to  installation  of  the  final  side  plate. 

The  wall  upon  which  the  primary  beam  was  incident,  and  all  walls  parallel  to 
it,  could  be  adjusted  in  thickness  from  1/2  in.  to  2.0  in. ,  allowing  a  span  of  areal 
density  range  normal  to  building  construction.  Several  measurements  of  the  steel 
used  gave  an  average  weight  of  41.2  psf  of  1  in.  steel  plate. 
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Figure  4.  Side  Plate  Being  Installed  on  Instrumented  Structure 


INSTRUME  NT  ATION 

The  majority  of  the  dose  measurements  were  made  with  uncompensated 
Victoreen  Model  362  (200  mr)  Pocket  Chambers,  positioned  by  insertion  with 
their  axes  perpendicular  to  the  primary  beam  in  thin  slabs  (1/4  in.)  of  Styrofoam 
of  a  density  of  0.033  g/cc.  A  preliminary  experiment  showed  that  the  response 
of  these  chambers  to  incident  ionizing  ladiation  was  essentially  flat,  provided 
that  the  chamber  axis  and  the  beam  direction  formed  an  angle  greater  than  about 
15°.  (For  smaller  angles,  the  sensitivity  dropped  only  about  15%. )  The  energy 
response  )f  this  instrument  has  been  shown  to  be  essentially  air-equivalent  down 
to  photon  energies  of  about  150  keV.^  The  decrease  in  charge  of  the  chamber 
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(which  is  proportional  to  exposure)  was  measured  with  a  Tech/Ops  Model  556 
Dosimeter  Charger-Reader,  a  portable  unit  designed  and  constructed  by  this  Com- 
pany  for  field  experiments.*^  The  charge  measurement  was  calibrated  against 
chamber  exposure  in  roentgens  by  use  of  a  standard  cobalt-60  source  on  a  calibra¬ 
tion  range.  Since  response  is  a  function  of  air  density,  corrections  were  made  for 
temperature  and  pressure. 

Where  radiation  levels  were  low,  readings  were  made  with  Victoreen  Model 
239A  Stray  Radiation  Chambers.  These  dosimeters,  operating  on  the  same  princi¬ 
ple  as  the  200  mr  pocket  chambers,  fully  discharge  with  an  accuracy  within  10% 
when  exposed  to  10  mr.  They  are  essentially  cylindrical— 4.5  in.  in  length  and 
2.0  in.  in  diam.  As  a  result,  where  radiation  levels  were  low,  total  mapping  of 
an  area  was  not  feasible  because  of  the  bulk  of  the  detector. 

A  rate  chamber  was  used  for  readings  requiring  extreme  sensitivity  but  no  pre¬ 
cise  positioning.  The  ionization  rate  chamber  used  has  a  volume  of  approximately 
320  cc  and  is  similar  to  the  one  found  in  a  standard  Victoreen  Model  592  Survey 
Meter.  The  chamber  signal,  which  is  provided  by  ion  currents  across  its  electrodes 
and  hence  is  proportional  to  the  exposure  rate,  enters  a  low-noise,  solid-state, 
feedback -type  amplifier  in  the  base  of  the  chamber.  The  amplified  output  is 
smoothed  through  a  simple  RC  circuit  and  fed  to  a  John  Fluke  digital  null  voltmeter. 
Exposure  rates  as  low  as  0.01  mr/hr  were  thus  read  with  excellent  reproducibility. 
The  rate  chamber  viras  required  only  for  readings  outside  the  egg  crate.  (Calibra¬ 
tion  of  this  instrument  also  requires  corrections  for  air  density,  through  the  param¬ 
eters  of  temperature  and  pressure.) 

Personnel  engaged  in  the  experiments  were  equipped  with  200  mr  direct-reading 
dosimeters  and  film  badges.  The  high  radiation  area  was  monitored  continuously 
v/ith  a  Tech/Ops  Model  492  Gammalarm  that  provided  a  backup  visual  indication  of 
the  source  position.  During  all  manipulation  of  the  structure  or  its  instrumentation, 
two  portable  survey  meters  for  personnel  monitoring  were  used. 

EXPERIMENTAL  PROCEDURE 

Before  any  experimental  data  reported  herein  were  taken,  the  cross-sectional 
outline  of  the  collimated  beam  was  determined. 
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First,  a  light  plywood  panel  was  installed  at  a  right  angle  to  the  source  pro¬ 
jector,  vertically  bisecting  the  projected  beam  at  a  distance  of  10  ft  from  the 
penumbral  collimator.  Charged  200-mr  pocket  chambers  were  attached  to  this 
target  in  both  vertical  and  horizontal  arrays.  The  source  was  "turned  on"  for  an 
8  min  exposure.  The  two  arrays  of  readings  were  used  to  provide  a  first  deter¬ 
mination  of  the  square  collimated  beam  position.  X-ray  film  was  then  attached  to 
the  plywood  target;  this  covered  the  region  around  the  outline  of  the  beam  as  deter¬ 
mined  by  the  chambers.  A  2  hr  exposure  was  made,  and  the  processed  film  gave 
a  clear  picture  of  the  beam.  To  determine  how  evenly  the  beam  intensity  was 
distributed  in  the  exposed  area  of  the  target,  film  badge  packets  were  arrayed  in  a 
cross  configuration  on  the  target  and  exposed  to  approximately  2  r.  The  configura¬ 
tion  and  readings  are  given  in  Figure  5. 


Figure  5.  Cross  Configuration,  with  Relative  Intensities  of 
Gamma  Radiation  Incident  on  Structure 

Upon  completion  of  the  preliminary  experiments  described  above,  the  box  ex¬ 
periments  were  run.  For  these  runs,  the  egg  crate  structure  was  placed  so  that 
for  normal  incidence  the  collimated  beam  just  covered  the  surface  of  the  front 
plate.  With  the  egg  crate  so  placed,  the  right  edge  of  the  front  plate  was  later 
used  as  the  pivot  axis  for  all  angular  experiments.  This  is  shown  in  Figure  6. 
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Figure  6 .  Compartmented  Structure  in  Beam  with  Pivot  Axis  to  Simulate 
Various  Angles  of  Incidence 


EXPERIMENTAL  RUN  AND  DATA  REDUCTION 

Each  run  involved  some  difference  in  structure  composition  or  orientation. 
Table  1  lists  these  briefly;  they  are  described  in  more  detail  here. 


TABLE  1 

CONSTRUCTION  OF  COMPARTMENTED  STRUCTURE 
FOR  EACH  EXPERIMENTAL  RUN 


Exp. 

No. 

Angle  of 
Incidence 

(<*eg) 

Front,  Back, 
and  Internal 

Side  Wall  Thickness 
(psf) 

Wall  Thickness 
(psf) 

Top 

Bottom 

Left 

Right 

1 

0 

20 

0 

0 

0 

0 

2 

0 

40 

0 

0 

0 

0 

3 

0 

80 

0 

0 

0 

0 

4 

0 

40 

40 

0 

0 

0 

5 

0 

40 

40 

40 

0 

0 

6 

0 

40 

40 

40 

40 

40 

7 

15 

40 

40 

40 

40 

40 

8 

30 

40 

40 

40 

40 

40 

9 

45 

40 

40 

40 

40 

40 

10 

60 

40 

40 

40 

40 

40 

11 

75 

40 

40 

40 

40 

40 

12 

0 

40 

40 

40 

0 

0 

13 

0 

40 

40 

40 

0 

40 

14 

45 

40 

40 

40 

40 

40 

15 

0 

40 

0 

0 

0 

0 

16 

0 

1  « 

0 

0 

0 

0 

Experiments  1  through  3  were  performed  with  no  side,  top,  and  bottom  plates. 
Thus  the  configuration  is  that  of  four  2-ft  square  steel  plates  spaced  1  ft  apart . 

The  incident  beam  exactly  covered  the  first  plate.  The  thickness  of  the  plates  was 
1/2  in.  in  the  first  three  experiments,  and  it  was  doubled  for  each  subsequent  series 

of  experiments . 
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Experiments  4  through  6,  all  performed  with  1  in.  plates,  were  a  series  de- 
signr'  1  to  investigate  the  effects  on  the  radiation  field  of  adding  plates  that  closed 
in  the  compartments.  The  top  plate  was  added  first,  then  the  bottom,  and  then 
the  sides. 

Experiments  7  through  1 1  were  performed  on  the  fully  completed  structure  and 
were  designed  to  investigate  the  effects  of  slant  incidence  of  the  beam.  Figure  7 
shows  the  compass  rose  drawn  on  the  floor  for  correct  positioning  of  the  structure, 
so  that  the  projector  beam  would  strike  the  front  plate  at  the  correct  angle  but 
would  not  strike  the  side  plates. 


Figure  7 .  Compass  Rose  Constructed  on  Floor  to  Facilitate 
Positioning  for  Angular-Incidence  Runs 
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Experiments  12  and  13  were  run  to  enable  us  to  estimate  the  dose  scattered 
into  a  compartment  from  the  side  walls  with  the  beam  at  normal  incidence  to  the 
structure.  For  this  series  of  experiments  we  arrayed  twenty-three  chambers  at 
1  in.  intervals  so  that  they  horizontally  traversed  the  center  of  the  front  com¬ 
partment.  The  chambers  were  positioned  vertically.  Experiment  12  was  run 
without  side  plates,  while  Experiment  IS  was  run  with  one  side  plate  installed. 

Thus ,  the  difference  in  dose  rates  measured  at  each  position  for  these  two  experi¬ 
ments  can  be  attributed  to  scatter  from  the  side  plate. 

Experiment  14  was  run  to  check  the  accuracy  with  which  we  could  place  the 
structure  for  runs  at  other  than  normal  incidence  to  the  beam.  The  same  array 
of  chambers  was  used  as  in  Experiments  12  and  13.  The  box  structure  was  placed 
at  45°  incident  to  the  beam. 

Experiment  15  was  run  to  ascertain  whether  any  detectable  stray  radiation  at 
the  outer  fringe  of  the  collimated  beam  could  scatter  into  the  front  compartment 
upon  striking  high  density  material.  Specifically,  this  would  indicate  if  a  correc¬ 
tion  would  have  to  be  made  for  scattering  off  the  legs  of  the  box-type  structure . 

This  experiment  was  essentially  a  repeat  of  Experiment  2 ,  but  a  steel  bar  (with  a 
1  in.  square  cross  section)  was  placed  on  top  of  the  front  plate.  Any  significant 
deviation  between  the  results  of  this  experiment  and  Experiment  2  would  indicate 
that  the  beam  overlapped  the  target,  and  scattering  of  the  fringe  radiation  would 
be  significant .  Plotting  the  dose  rates  against  the  position  of  each  chamber  should 
result  in  a  curve  that  breaks  steeply  at  the  edge  of  the  direct  besun.  This  would 
give  a  good  indication  of  the  accuracy  of  alignment  of  the  box-type  structure  with 
the  collimated  beam . 

Experiment  16  was  a  phantom  structure  experiment  run  primarily  as  a  cali¬ 
bration.  The  box  structure  was  removed,  and  in  its  place  a  light  plywood  panel 
was  positioned  to  correspond  with  the  front  plate  of  the  box.  We  then  cut  out  of 
the  panel  a  2  sq  ft  opening  aligned  in.  the  path  of  the  collimated  beam .  Masking 
tape  was  used  to  install  nine  pocket  chambers  in  this  opening,  at  positions  corre¬ 
sponding  to  those  used  during  box  experiments  and  at  a  depth  corresponding  to  the 
front  surface  of  the  front  plate  of  the  box  structure  (at  the  normal  angle  of  incidence) . 
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Experiments  17  and  18  were  a  further  attempt  to  read  scatter  back  into  the 
box  from  a  side  plate  with  the  beam  at  a  normal  incidence .  For  this  series ,  the 
ionization  rate  chamber  was  placed  at  points  1  through  5  as  shown  in  Figure  8. 
Readings  were  taken  with  both  sides  removed  from  the  box  and  also  with  alternate 
sides  installed. 


BEAM 


REMOVE  FOR  EXP.  18 


FRONT 


CENTER 


REAR 


O  O  O  O  O 

I  2  3  4  S 

Figure  8.  Positioning  of  Rate  Chamber  for  Experiments  17  and  18 


All  data  given  in  the  tables  are  in  terms  of  r/hr,  assuming  a  1  r/hr  dose  rate 
at  the  position  with  the  structure  removed.  Thus,  all  data  may  be  reduced  to  ratios 
of  exposure  rates  of  interest  to  the  exposure  rate  in  the  incident  beam.  Chamber 
readings  were  converted  to  dose  rates  by  multiplying  by  a  constant  that  was  deter¬ 
mined  by  calibration  against  a  standard  source  in  air.  The  calibration  constant, 
plus  the  temperature  and  pressure  corrections  postulated  by  air-wall  ionization 
chamber  theory,  relate  the  chamber  e2q;x}sure  to  the  readout  measurement  /za  by 
the  formula  . 


D  = 


8.65P 


where  T  is  in  degrees  Rankine,  P  in  inches  of  mercury,  and  D  in  milliroentgens. 
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EXPERIMENTAL  DATA 

Unless  otherwise  indicated,  ion  chambers  were  positioned  in  each  of  the 
three  egg  crate  compartments  in  three  vertical  planes,  with  nine  chambers  to  a 
plane .  Thus  this  total  of  eighty-one  readings  per  exposure  provided  excellent  de¬ 
tail  of  the  radiation  field.  The  compartments  are  identified  as  front,  center,  and 
rear,  with  reference  to  the  nearness  to  the  radiation  source;  the  instrument  planes 
are  similarly  identified  as  2,  6,  and  10,  representing  their  distances  in  inches  from 
the  front  of  the  respective  compartment.  Detectors  in  each  plane  were  spaced  8  in. 
apart;  these  are  identified  by  the  letter  code  shown  in  Figure  9.  The  planes  are 
sketched  in  Figure  10.  Note  that  detectors  E,  O,  and  Y  are  closest  to  the  pivot 
axis  used  in  Experiments  11  and  14. 


Dividing  measured  chamber  ex¬ 
posure  by  exposure  time  results  in 
measured  dose  rates.  These  were  cor¬ 
rected  for  temperature  and  pressure 
effects  as  mentioned  above.  The  dose 
rates  thus  obtained  were  divided  by  dose 
rates  measured  at  the  location  of  the 
front  slab  of  the  structure  (with  the 
structure  absent).  This  ratio  is,  of 
course,  affected  by  the  source-detector 
distance,  and  to  correct  for  this,  a  fac¬ 
tor  of  the  ratio  of  the  squares  of  the  dis¬ 
tances  of  the  detector  and  the  front  plate  was  used.  Finally,  a  1.05  factor  was  used 
to  correct  for  the  presence  of  about  5%  soft  radiation  (backscatter  or  X-rays  from 
the  projector),  which  was  suggested  by  later  single-slab  measurements.  With  these 
corrections,  the  basic  data  are  presumed  to  be  normalized  to  a  monoenergetic 
parallel-ray  beam  of  radiation  that  has  unit  intensity  as  it  strikes  the  structure 
being  examined. 

The  basic  data  obtained  in  Experiments  1  through  11  and  corrected  as  stated 
above  are  presented  in  Tables  2  through  12. 
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Figure  9.  Detector  Plane  Showing 
Letter  Code 


TABLE  3 


EXPERIMENT  2:  ATTENUATION  BY  40  PSF  WALLS, 
WITHOUT  TOP,  SIDES,  AND  BOTTOM 


Compart¬ 

ment 

Detector 

Plane 

(in.) 

Detector  Position 

A 

c 

E 

K 

M 

O 

U 

W 

V 

2 

0.55 

0.56 

0.54 

0.59 

0.59 

0.59 

0.58 

0.58 

0.60 

Front 

6 

.50 

.52 

.49 

.52 

.54 

.52 

.49 

.50 

.50 

10 

.45 

.48 

.45 

.48 

.51 

.48 

.46 

.48 

.46 

2 

.24 

.24 

.24 

.25 

.25 

.25 

.24 

.24 

.23 

Center 

6 

.21 

.23 

.21 

.23 

.24 

.23 

.21 

.23 

.19 

10 

.19 

.21 

.19 

.21 

.22 

.21 

.19 

.22 

.18 

2 

B 

.097 

.094 

.109 

.097 

.096 

.094 

.092 

Rear 

6 

.093 

.084 

.088 

.089 

.081 

.083 

10 

.087 

.078 

.082 

.088 

.082 

.076 

.084 

.075 

TABLE  4 

EXPERIMENT  3:  ATTENUATION  BY  80  PSF  WALLS, 
WITHOUT  TOP,  SIDES,  AND  BOTTOM 


Compart- 

Detector 

Plane 

(in.) 

Detector  Position 

ment 

A 

C 

E 

K 

M 

o 

u 

w 

V 

2 

0.30 

0.31 

0.30 

0.31 

0.32 

0.31 

0.31 

0.31 

0.31 

Front 

6 

.26 

.28 

.26 

00 

• 

.31 

.30 

.26 

.27 

.26 

10 

.22 

.27 

.26 

.26 

.29 

.26 

.23 

.26 

.23 

2 

.052 

.057 

.054 

.057 

.060 

.056 

.052 

.054 

.052 

Center 

6 

.046 

.052 

.048 

.049 

.056 

.050 

.045 

.050 

.045 

10 

.041 

.047 

.042 

.044 

.051 

.046 

.041 

.046 

.042 

2 

.0097 

.0104 

Rear 

6 

.0087 

10 

.0083 

.0089 

18  B  U  K  L 

1  N  0 

T  0 

N  • 

M 

A  t  S 

A  C  H 

u  s 

t  T  T 

3 

TABLE  5 


EXPERIMENT  4:  ATTENUATION  BY  40  PSF  WALLS,  WITH  40  PSF  BOTTOM 


Compart- 

Detector 

Plane 

(in.) 

Detector  Position 

meat 

A 

C 

E 

K 

M 

O 

U 

W 

V 

2 

0.55 

0.58 

0.55 

0.61 

0.61 

0.61 

0.59 

0.61 

0.60 

Front 

6 

.50 

.52 

.50 

.53 

.56 

.53 

.52 

.54 

.50 

10 

.44 

.49 

.46 

.50 

.52 

.48 

.48 

.50 

.47 

2 

.24 

.24 

.25 

.25 

.27 

.25 

.24 

.25 

.23 

Center 

6 

.21 

.21 

.21 

.21 

.25 

.23 

.20 

.24 

.20 

10 

.19 

.21 

.21 

.21 

.22 

.21 

.21 

.22 

.19 

2 

.099 

.108 

.099 

.101 

.116 

.106 

.098 

.109 

.102 

Rear 

6 

.085 

.096 

.085 

.092 

.103 

.094 

.081 

.098 

.094 

1 

10 

.083 

.087 

.079 

.083 

.093 

.083 

.075 

.089 

.083 

TABLE  6 

EXPERIMENT  5:  ATTENUATION  BY  40  PSF  WALLS, 
WITH  40  PSF  TOP  AND  BOTTOM 


Compart- 

Detector 

Plane 

(in.) 

Detector  Position 

ment 

A 

C 

E 

K 

M 

o 

u 

w 

V 

2 

0.54 

0.58 

0.56 

0.59 

0.61 

0.60 

0.58 

0.59 

0.58 

Front 

6 

.49 

.52 

.48 

.52 

.55 

.53 

.49 

.52 

.50 

10 

.45 

.49 

.45 

.48 

.51 

.48 

.45 

.49 

.46 

2 

.23 

.25 

.24 

.24 

.25 

.24 

.23 

.23 

.22 

Center 

6 

.21 

.23 

.21 

.20 

.24 

.21 

.20 

.21 

.20 

10 

.19 

.22 

.19 

.19 

.23 

.21 

.18 

.21 

.18 

2 

.096 

.102 

.094 

.096 

.107 

.099 

.094 

.103 

.096 

Rear 

6 

.079 

.093 

.083 

.087 

.094 

.087 

.079 

.093 

.087 

10 

.080 

.088 

.076 

.080 

.087 

.080 

.072 

.087 

.076 
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TABLE  7 


EXPERIMENT  6:  ATTENUATION  BY  40  PSF  COMPLETE  STRUCTURE ,  WITH 
40  PSF  TOP,  BOTTOM,  SIDES,  AND  WALLS 


Cc^nijart- 

Detector 

Plane 

(in.) 

Detector  Position 

ment 

A 

C 

E 

K 

M 

O 

U 

W 

V 

2 

0.55 

0.58 

0.55 

0.59 

0.59 

0.59 

0.58 

0.59 

0.59 

Front 

6 

.50 

.53 

.49 

.53 

.54 

.55 

.50 

.53 

.50 

10 

.46 

.49 

.47 

.49 

.51 

.48 

.48 

.49 

.47 

2 

.23 

.25 

.25 

.25 

.27 

.25 

.23 

.23 

.23 

Center 

6 

.21 

.24 

.23 

.21 

.24 

.23 

.20 

.21 

.20 

10 

.21 

.22 

.21 

.21 

.23 

.21 

.19 

.21 

.19 

2 

.098 

.102 

.098 

.098 

.112 

.102 

.096 

.106 

.098 

Rear 

6 

.087 

.093 

.087 

.089 

.098 

.090 

.080 

.097 

.089 

10 

.080 

.088 

.078 

.082 

.091 

.082 

.074 

.091 

.078 

TABLE  8 

EXPERIMENT  7:  ATTENUATION  BY  40  PSF  COMPLETE  STRUCTURE, 
WITH  BEAM  INCIDENT  AT  15°  FROM  NORMAL 


Compart- 

Detector 

Plane 

(in.) 

Detector  Position 

ment 

A 

C 

E 

K 

M 

0 

u 

w 

V 

2 

0.55 

0.55 

0.53 

0.59 

0.58 

0.56 

0.55 

0.56 

0.55 

Front 

6 

.54 

.53 

.45 

.57 

.54 

.47 

.50 

.50 

.47 

10 

.50 

.50 

.32 

.52 

.52 

.35 

.48 

.48 

.31 

2 

.25 

.25 

.13 

.27 

.25 

.14 

.24 

.23 

.12 

Center 

6 

.24 

.23 

.087 

.24 

.24 

.094 

.21 

.21 

.076 

10 

.22 

.22 

.062 

.22 

.22 

.062 

.21 

.19 

.053 

2 

.098 

.096 

.023 

.109 

.106 

.025 

.099 

.096 

.021 

Rear 

6 

.089 

.085 

.021 

.096 

.090 

.024 

.085 

.084 

.021 

10 

.084 

.075 

.018 

.089 

.082 

.022 

.079 

.075 

.018 

20  ■ 

tl  K  L 

1  N  0 

T  0 

N  • 

M  A 

s  t 

A  C  H 

u  s 
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EXPERIMENT  8:  ATTENUATION  BY  40  PSF  COMPLETE  STRUCTURE , 
WITH  BEAM  INCIDENT  AT  30°  FROM  NORMAL 
(See  Figure  7) 


Compart- 

Detector 

Plane 

(in.) 

Detector  Position 

ment 

C 

E 

K 

M 

O 

u 

W 

V 

2 

0.50 

0.52 

0.39 

0.52 

0.53 

0.40 

0.50 

0.52 

0.40 

Front 

6 

00 

.47 

.15 

.46 

.48 

.19 

.47 

.46 

.16 

10 

.44 

.44 

.072 

.48 

.47 

.078 

.44 

.42 

.065 

2 

.19 

.021 

.21 

.20 

.024 

.20 

.17 

.020 

Center 

6 

.18 

.114 

.021 

.19 

.12 

.025 

.17 

1C 

.16 

.055 

.018 

.17 

.068 

.023 

.16 

.052 

.021 

2 

.069 

.020 

Rear 

6 

.018 

10 

.034 

.013 

.0075 

.038 

.014 

.0085 

■SI 

.013 

TABLE  10 

EXPERIMENT  9;  ATTENUATION  BY  40  PSF  COMPLETE  STRUCTURE, 
WITH  BEAM  INCIDENT  AT  450  FROM  NORMAL 


Compart¬ 

ment 

Detector 

Plane 

(in.) 

Detector  Position 

A 

c 

E 

K 

M 

0 

u 

W 

V 

2 

0.39 

0.42 

0.32 

0.42 

0.42 

0.33 

0.39 

0.40 

0.33 

Front 

6 

.38 

.36 

.069 

.37 

.38 

.076 

.36 

.36 

.060 

10 

.35 

.26 

.045 

.37 

.26 

.049 

.35 

.23 

.040 

2 

- 

— 

.067 

— 

— 

.064 

— 

— 

Center 

6 

.025 

.013 

.049 

.027 

.014 

.046 

.023 

.0113 

10 

- 

- 

.021 

- 

— 

.023 

— 

— 

Rear 

6 

.0087 

.0065 

.0048 

.0099 

.0071 

.0052 

.0087 

.0063 

.0048 
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TABLE  11 

EXPERIMENT  10:  ATTENUATION  BY 
40  PSF  COMPLETE  STRUCTURE,  WITH 
BEAM  INCIDENT  AT  60°  FROM  NORMAL 


Compart¬ 

ment 

Detector 

Plane 

(in.) 

Detector  Postion 

A 

C 

E 

K 

M 

O 

U 

W 

V 

2 

0.23 

0.23 

0.034 

0.24 

0.26 

0.037 

0.23 

0.26 

0.034 

Front 

6 

.21 

.084 

.021 

.21 

.099 

.024 

.21 

.084 

.021 

10 

.101 

.042 

.018 

.12 

.042 

.021 

.093 

.042 

.018 

Center 

6 

.015 

.0106 

.0106 

.018 

.013 

.0111 

.014 

.0102 

.0102 

Rear 

6 

.0034 

.0024 

.0027 

.0037 

.0030 

.0030 

.0032 

.0024 

.0026 

TABLE  12 

EXPERIMENT  11:  ATTENUATION  BY 
40  PSF  COMPLETE  STRUCTURE,  WITH 
BEAM  INCIDENT  AT  75°  FROM  NORMAL 


Compart¬ 

ment 

Detector 

Plane 

(in.) 

Detector  Position 

a 

c 

E 

K 

M 

O 

u 

w 

V 

2 

["■  . . . 

0.063 

- 

— 

0.069 

- 

-  0 

.064 

- 

— 

Front 

6 

.034 

0.021 

0.0058 

.036 

0.023 

0.0067 

.033 

0.019 

0.0055 

10 

.020 

- 

— 

- 

- 

- 

.020 

- 

— 

Center 

6 

.0038 

.0055 

.0039 

.0043 

.0061 

.0052 

.0036 

.0052 

.0046 

Rear 

6 

.00097 

.0014 

.0014 

.0107 

.0014 

.0014 

.00101 

.0014 

.0014 
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Figure  11.  Horizontal  Traverse 
of  Radiation  Field  in  Midplane  of 
Each  Compartment  (Experi¬ 
ments  12  and  13) 


Results  of  Experiments  12  and 
13,  designed  for  a  more  detailed 
investigation  of  the  field  profile  in 
the  structure  normal  to  the  beam, 
are  graphed  in  Figure  11 .  Since  the 
chambers  were  only  1  in.  apart,  it 
was  possible  to  see  the  edge  ejects 
as  the  side  of  the  structure  was 
approached.  Incidentally,  this  fig¬ 
ure  provides  a  demonstration  of  the 
accuracy  of  the  instrumentation  in 
that  no  points  deviate  more  than 
about  4%  or  5%  of  their  positions 
on  a  smooth  curve  joining  their 
neighbors . 


Experiments  14  and  15 -essen¬ 
tially  checks  on  the  accuracy  of 
positioning  of  the  structure -gave 
negative  results;  that  is,  no  detec¬ 
table  differences  existed  in  these 
configurations  from  dose  rates  re¬ 
ceived  in  similar  experiments 
run  previously .  (Fundamentally , 
Experiments  14  and  15  were  repeats 
of  Experiments  12  and  2  respectively. 


Finally,  Experiment  16  showed  that  the  front  of  the  structure  received  a  dose 
rate  of  1.09  r/hr.  This  figure  was  used  to  determine  the  ratios  reported  in  Tables 
2  through  12. 
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CHAPTER  4 


ANALYSIS  OF  DATA 


ATTENUATION  BY  FINITE  SLABS  OF  NORMAJ.LY  INCIDENT  RADIATION 


The  gamma  energy  flux  resulting  from  photons  that  have  undergone  no  inter¬ 
action  since  their  birth  and  impinging  on  a  small  isotropic  detector  is  given  by 
the  integral  over  all  space  of  the  source  density  (energy /sec /cc)  at  a  distance  r, 
divided  by  the  quantity 


2 

4irr  exp 


r 


M(r)dr  , 


(1) 


where  /i(r)  is  the  linear  attenuation  coefficient  of  the  medium  at  the  point  being 
traversed  by  the  photon.  The  attenuation  coefficient  depends  on  the  photon 
energy  and  the  density  and  nature  of  the  substances  being  traversed . 

Rarely  does  this  calculated  value  of  the  energy  flux  due  to  uncollided  photons 
approach  the  true  flux  at  the  detector .  Scattered  radiation  usually  increases  the 
response  of  the  detector  by  a  factor  B,  referred  to  as  the  "buildup  factor"  for  the 
effect  in  question.  Exact  values  of  B,  which  are  generally  for  simple  geometries 
such  as  point  sources  in  infinite  media,  are  obtained  with  difficulty  only  by  ardu¬ 
ous  machine  calculations . 

As  an  example,  consider  that  the  centerline  detector  responses  in  the  egg 
crate  have  normal-beam  incidence.  This  apparently  simple  geometry— four 
square  slabs  with  a  beam  of  parallel  radiation  striking  the  outer  surface  of  the 
first-provides  experimental  values  of  B  that  are  to  be  analyzed  in  the  light  of 
previous  calculations. 

Placing  a  detector  in  the  center  of  the  first  compartment,  we  get  a  response 
(Table  3)  that  is  0.54  of  that  obtained  with  the  egg  crate  removed .  The  1  in .  slab 
provides  1.07  mfp  to  the  ray;  this  results  in  a  narrow  beam  attenuation  factor  of 
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0.343.  Thus,  the  value  of  B  is  1.525;  the  additional  52.5%  must  be  attributed  to 
radiation  that  either 

(1)  penetrates  the  slab  as  scattered  radiation,  or 

(2)  rebounds  from  the  posterior  slabs  as  an  albedo  effect. 

A 

The  energy  albedo  of  iron  for  cobalt  radiation  is  about  0.025.  The  albedo 
depends,  of  course,  on  the  angular  distribution  of  the  incident  radiation  and  on  the 
thickness  of  the  iron,  but  the  precise  size  of  the  effect  is  obviously  of  little  im¬ 
portance  and  may  be  taken  into  account  here  by  accepting  this  figure  as  the  albedo 
0! .  Most  of  the  observed  buildup  is  therefore  due  to  process  (1) . 

Spencer  has  reported  a  number  of  calculations  of  the  attenuation  of  cobalt 
radiation  through  concrete.  He  defines  a  function  spc,  6^)  that  gives  the  dose 
attenuation  by  concrete  of  X  psf  thickness  for  radiation  from  a  source  plane  in¬ 
finite  in  extent,  with  all  radiation  leaving  the  plane  at  an  angle  9^  to  the  normal. 
Now  s(X,  6^)  readily  provides  us  with  an  upper  limit  to  the  detector  response,  to 
be  achieved  only  u4ien  the  iron  slabs  are  infinite  planes,  while  the  exponential  func¬ 
tion  provides  us  with  a  lower  limit  that  is  approached  only  as  the  lateral  extent  of 
slabs  and  detector  vanishes  to  zero.  Table  13  gives  these  limits,  together  with 
centerline  detector  responses,  and  shows  that,  indeed,  the  pragmatic  case  of  finite 
slabs  does  fall  between  the  limits. 


TABLE  13 

COMPARISON  OF  NARROW-BEAM,  WIDE-BEAM,  AND 
FINITE -BEAM  ATTENUATION  BY  IRON 


Mt 

Narrow- Be  am 

Q-Pt 

Wide- Beam 
s(X) 

Measured 

D(X) 

1.07 

0.343 

0.690 

0.54 

2.14 

.118 

.385 

.24 

3.21 

.040 

.160 

.097 

i  u  II  1. 
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Since  the  spread  in  predictions  is  much  too  wide,  an  expression  for  the 
buildup  that  incorporates  a  correction  for  the  effect  of  slab  extent  is  required. 

One  method  for  handling  this  problem  is  the  Monte  Carlo  calculation^  to  be 
discussed  in  Chapter  5 .  The  results  compare  fairly  closely  with  the  experimental 
values.  For  a  slab  thickness  of  1.07  mfp,  the  result  calculated  by  the  Monte  Carlo 
method  is  0.558 -a  difference  from  experimental  values  of  less  than  2%. 

ATTENUATION  BY  FINITE  SLABS  OF  ISOTROPIC  RADIATION 

From  the  results  given  for  the  parallel-beam  case  above,  we  expect  similar 
behavior  trends  for  the  case  where  the  incident  radiation  is  isotropic.  The  trends 
are  more  complex,  since  now  various  angles  of  incidence  are  involved.  One 
would  expect  that  the  extent  of  the  slab  would  still  affect  the  measured  dose  rate, 
however,  and  that  a  construction  of  the  buildup  function  from  angular  and  depth 
components  would  be  possible .  The  Engineering  Manual  presents  theoretical 
comparisons  for  the  results  of  such  a  construction.  The  carves  presented  therein 
for  the  evaluation  of  rooftop  contamination  are  formed  by  the  product  L(X)L  pC,  co) 

V 

obtained  by  Spencer.  In  such  a  case,  L(X)  is  the  attenuation  of  the  isotropic  in¬ 
finite-layer  source  by  material  of  thickness  X  and  corresponds  to  a  more  precise 
evaluation  of  the  dose  rate  (this  is  for  an  infinite  slab  and  an  isotropic  source'.), 
while  L  (X,  cj)  is  the  fraction  of  the  radiation  at  X  resulting  from  the  part  of  the 
plane  source  that  subtends  the  solid  angle  fraction  w  at  the  detector. 

The  methods  presented  do  not,  however,  make  provision  for  the  air  spaces 
between  detector  and  slab  and  for  the  division  of  the  slab  into  layers.  Therefore, 
it  is  of  interest  to  compare  results  obtained  experimentally  with  results  obtained 
from  curves  of  LL  .  Although  this  cannot  be  done  for  the  parallel-beam  situation, 
it  can  be  done  for  the  isotropic  beam. 

To  simulate  an  isotropic  beam,  we  directed  the  parallel  beam  at  the  egg 
crate  at  the  different  angles  '/ith  the  normal,  as  discussed  in  Chapter  2.  It  was 
then  possible  to  perform  a  summation,  suitably  weighted  over  the  angles  of  inci¬ 
dence,  of  the  dosimeter  exposures  obtained  and  to  compare  this  with  predictions  for 
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isotropic  emission  from  the  face  of  the  first  slab.  Because  this  summation 
had  to  be  done  with  care,  the  method  described  below  was  derived  for  a  correct 
weighting  and  summing  of  the  exposures . 

CONSTRUCTION  OF  ISOTROPIC  SOURCE  RESPONSE  FROM 
MEASUREMENTS  OF  PARALLEL  BEAM  RESPONSE 

For  each  angle  O  that  the  rays  make  with  the  normal  to  the  slabs ,  there  is  a 
range  of  azimuthal  angles  <f>  that  contributes  no  direct  beam  to  a  given  detector . 
For  a  simple  example  of  this  forbidden  range,  consider  a  detector  in  the  center 
of  the  middle  compartment.  From  Figure  12,  it  can  be  seen  that  rays  with  angles 
0  less  than  or  equal  to  0.^  strike  the  detector  with  all  angles  <f>,  but  rays  with 
angles  greater  than  this  are  restricted  to  four  ranges  of  <p  that  correspond  to  the 
corners  of  the  top  plate.  Finally,  no  rays  of  angles  0  equal  to  or  greater  than 
may  strike  this  detector  without  having  first  scattered. 


SOURCE  SLAB 


Figure  12.  Allowed  Angles  of  0  for  Straight-Line 
Propagation  of  Radiation 
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Figure  13.  Construction  for  Deter¬ 
mination  of  Eq.  3 


The  fraction  of  acceptable  angles  <t> 
to  the  total  range  of  angles  can  be  deter¬ 
mined  by  geometric  considerations  for 
the  cases  between  6^  and  6^  for  center- 
line  detectors.  Consider  the  circle  of 
loci  of  sources  positioned  on  the  front 
plate  to  emit  at  angle  6  (and,  arbitrary, 
<^)  so  as  to  strike  the  detector  with  direct 
rays .  The  fraction  sought  is  then  the 
fraction  of  the  circle  that  is  on  the  plate . 
With  the  notation  of  Figure  13,  this  is 
equal  to  (4L)/(27R);  it  can  also  be  seen 
from  the  figure  that 


sin  4  =  ^  , 

(L+M)  =  4R, 


and 


L+  2M  =  I  R. 


With  these  identities,  the  range  of  <p  is  given  by 


2ir 


p(0 , 0)d<p  =  ^  arc  sin  cot  0^-1, 


(t  - 


(2a) 

(2b) 


(2c) 


(3) 


where  2A  is  the  length  of  a  side  of  the  square  plate,  and  Z  is  the  distance  of  the 
detector  from  the  plate.  The  range  is  plotted  in  Figure  14  for  detectors  in  the 
centers  of  each  of  the  three  compartments. 
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Figure  14.  Fraction  of  the  Azimuthal  Angle  Range  That  Contributes 
Direct-Beam  Radiation  to  a  Detector  in  the  Center  of  a 
Given  Compartment 

As  stated  in  the  beginning  of  this  chapter,  the  direct-beam  component  of  the 
radiation  is  attenuated  exponentially,  while  the  scattered  component  behaves 
according  to  much  more  complex  laws .  Therefore ,  we  assume  our  knowledge  of 
the  direct-beam  component  is  correct;  in  particular,  we  subtract  from  the  experi¬ 
mental  data  an  analytic  function  corresponding  to  the  direct  beam.  This  function 
has  a  sharp  cutoff,  as  discussed  above.  The  residual  component,  however,  can  be 
expected  to  behave  fairly  smoothly,  which  allows  us  to  integrate  it  numerically 
over  all  angles  using  a  small  number  of  intervals .  (Because  the  detectors  are  not 
at  ideal  points ,  the  cutoff  is  often  not  sharp  enough;  thus  some  angles  provide 
measurements  at  some  detectors  that  are  not  usable  in  this  analysis . ) 

An  analytic  function  describing  the  direct  beam  at  a  centerline  detector  due  to 
radiation  incident  on  the  front  plate  with  angle  9  to  the  normal  is  exp  jj-  X  sec  o]  , 
where  X  is  the  mean-free-path  thickness  of  slabs  in  front  of  the  detector .  Integration 
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of  this  function  over  a  source  area  gives  the  unscattered  beam  dose  rate  D  for 

u 

an  isotropic  rooftop  source  as 


D  = 
u 


r  dA  ^ 

A 


(4) 


For  a  detector  at  a  distance  z  from  a  square,  with  a  shielded  source  area  of 
size  2a,  the  unscattered  component  is 


2^2 
r  +  z 


a 


(5) 


where  the  second  integral  represents  the  component  due  to  the  portion  of  the  source 
plane  outside  a  concentric  circle  of  diameter  2a,  and  P(6)  is  the  function  graphed 
in  Figure  14.  The  first  component  of  is  readily  integrable  in  terms  of  tabulated 
functions  and  is 


E^(X)  - 


(6) 


The  second  integral  must  be  evaluated  numerically . 

To  obtain  the  total  dose  rate  at  the  detector  positions ,  we  must  add  a  scattered- 
beam  dose  component  to  the  values  of  obtained  above.  This  may  be  accom¬ 
plished  by  defining  D  (cos  9)  as  the  scattered-beam  dose  component  due  to  source 
radiation  at  an  angle  9  with  the  centerline  of  the  structure.  Then  the  total  dose  is 


D  =  D  +  D  , 
u  s  ’ 


(7) 
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and  the  scattered-beam  dose  component  is 


1 

Dg  =  y  Dg(cos  0)  d(cos  0)  . 
o 


(8) 


Values  of  Dg(cos  0)  are  obtained  experimentally.  We  subtract  the  unscattered 

mm  SGC  ^ 

component  e  from  the  total  measured  dose  rates  of  Tables  7  through  12  for 

centerline  positions.  The  resultant  values  must  be  divided  by  the  cosine  of  the  in¬ 
cidence  angle,  since  the  intensity  of  radiation  (current  per  unit  area  of  front  plate) 
was  experimentally  proportional  to  cos  0 .  The  results  of  this  manipulation  gener¬ 
ate  the  family  of  curves  given  in  Figure  15.  The  integrals  under  these  curves  are 

,  and  are  recorded  in  Table  14  for  the  nine  centerline  detectors . 
s 


cos  9 


Figure  15.  Scattered  Component  of  Dose  Rates  Along  Centerline  of 
Compartments,  Due  to  Incident  Radiation  and  Angles  0 
with  Normal  to  Front  Plate 
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TABLE  14 

COMPONENTS  OF  ISOTROPIC  ROOF  SOURCE  DOSE  CONTRIBUTION 


Compartment 

Detector 

Plane 

(in.) 

D 

u 

D 

s 

4i 

D(experimental) 

2 

0.2004 

0.2846 

0.4850 

Front 

6 

.1714 

.1807 

.3521 

10 

.1246 

.1047 

.2293 

2 

.0276 

.0736 

.1012 

Center 

6 

.0204 

.0423 

.0627 

10 

.0196 

.0277 

.0473 

2 

.0040 

.0257 

.0297 

Rear 

6 

.0034 

.0129 

.0163 

10 

.0028 

.0061 

.0089 

Ratio  of  dose  rate  with  40  psf  floors  to  dose  rate  with  0  psf  floors . 


COMPARISON  WITH  ROOFTOP  CONTAMINATION  CALCULATIONS 

2 

The  OCD  Engineering  Manual  presents  a  method  for  the  calculation  of  the 
effects  of  radioactive  contamination  on  a  rooftop.  The  results  obtained  in  the 
previous  chapter  and  tabulated  above  in  Table  14  should  be  comparable  to  those 
obtained  by  this  method. 

We  calculate  that  a  1  in.  steel  plate  provides  an  attenuation  equivalent  to 
38.4  psf  of  water.  From  this,  and  from  a  calculation  of  the  solid  angle  sub¬ 
tended  by  the  roof  at  the  detector,  we  obtain  the  parameters  X  (equivalent  pounds 
per  square  foot)  and  co  (subtended  solid-angle  fraction) .  The  dose  rate  due  to 
roof  sources  is  then  D(a; ,  X) . 
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With  vanishing  X,  the  dose  rate  due  to  roof  sources  is  readily  calculated 
from  an  integral.  Where  the  source  concentration  is  such  that  current  from  a 
unit  area ,  in  directions  about  a  solid  angle  dO ,  gives  a  dose  rate  equal  to 
(di2/27r),  the  expression  for  the  dose  rate  is 


plane 


dxdy 

2 

27rp 


(9) 


where  p  is  the  distance  from  the  unit  area  (x,  y)  to  the  detector.  This  integral  is 
readily  evaluated  if  the  plane  is  a  disk  and  the  detector  is  located  along  the  disk 
axis .  The  dose  is  then 


D(without  shield)  =  -£n(l  -  w)  . 


(10) 


To  obtain  the  reduction  factor  R^,  one  needs  to  know  the  ratio  of  the  dose  rate  due 
to  shielded  sources  in  all  directions,  to  the  dose  rate  due  to  the  same  sources 
when  they  are  unshielded .  The  shielded  dose  rate  may  be  decomposed  into  con¬ 
tributions  from  sources  in  various  solid  angle  fractions,  so  that 


«f  = 


D(aJ^,X^)  +  D(u;2,X2)  + 


(11) 


where  X^  is  the  shielding  provided  by  3  ft  of  air.  This  ratio  is  obviously  equal  to 


■D(a;^,X^)' 

A 

'  D(‘^2*^2)' 

D  (w, ,  X  ) 
o'  1  o' 

_Do(2.X^)_ 

•r 

(12) 


where  Dq(w^i  X^)  may  be  chosen  as  the  dose  rate  due  to  sources  in  the  range  of 
c^i,  with  negligible  shielding  X^  present.  Then,  the  first  bracketed  ratio  in 
Eq.  (1^  may  be  computed  by  dividing  experimental  values  from  Table  14  by  values 
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computed  from  Eq.  (10).  The  second  ratio  is  equivalent  to  L^(X^,  as  dis¬ 
cussed  earlier  in  this  chapter  under  isotropic  radiation,  with  the  further  obser¬ 
vation  that  L(X^)  is  unity. 

Thus ,  the  reduction  factor  R^.  is  the  sum  of  all  the  dose  rate  contributions  C . 
The  overhead  contribution  is  given  in  the  OCD  Engineering  Manual.  In  pre¬ 
senting  this  function,  Eisenhauer  and  FitzSimons,  the  manual  authors,  utilized 

7 

Spencer's  curves;  specifically,  C  equals  the  product  L(K)L  (X,  w). 

The  present  experimental  data  are  not  exactly  applicable  to  a  test  of  , 
since  that  function  was  derived  for  gamma  radiation  of  a  spectrum  corresponding 
to  early  fallout.  However,  one  may  compute  a  similar  function  from  curves  of  L 
and  L  in  Ref.  7,  especially  for  cobalt-60  radiation  and  a  concrete  medium.  At 
these  gamma  energies,  the  mass  attenuation  coefficient  is  remarkably  insensitive 
to  the  nature  of  the  medium  (see,  for  example,  the  graph  in  Ref.  9),  and  data  for 
concrete  may  be  applied  to  measurements  in  steel. 

From  this  line  of  reasoning ,  it  might  be  presumed  that  our  data  are  not  too 
precise  a  test  of  the  accuracy  of  the  Spencer-Eisenhauer- FitzSimons  rooftop  con¬ 
tamination  curves.  However,  we  see  in  Table  15  that  several  figures  extracted 
from  Table  14  compare  well  with  the  dose  rates  derived  from  this  source.  The 
theory  consistently  evaluates  the  dose  rate  as  more  than  the  experiment  predicts . 
Even  if  we  could  assign  all  the  error  to  theory,  it  would  not  deviate  more  than 
Eisenhauer  had  anticipated,  according  to  Ref.  10. 


TABLE  15 

COMPARISON  OF  EXPERIMENTALLY  DERIVED  AND 
THEORETICAL  ROOFTOP-SOURCE  DOSE  RATES 


X 

0) 

♦ 

R 

d/dJ 

/TheorvN 
i  Exper .  ) 

38.4 

0.60 

0.468 

0.352 

1.33 

76.8 

.20 

.325 

.281 

1.15 

115.2 

.09 

.176 

.171 

1.03 

Ratios  of  L(X)L^(X,  w)/L(o)L^(o,  w)  from  Ref.  2,  p.  126. 
^D  from  Table  14,  D^  from  Eq.  (10). 
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The  comparison  indicates,  however,  that  a  worthwhile  effort  would  be  that  of 
measuring  the  actual  rooftop  radioactive  deposit  protection  factor  afforded  by 
the  roof  and  overhead  floors  of  a  building .  The  present  experiment  indicates  some 
departure  between  the  prediction  and  the  measurement,  possibly  due  to  separation 
distance  effects  or  to  the  separation  of  the  layers  (floors)  of  the  shield. 
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CHAPTER  5 

COMPARISON  WITH  MONTE  CARLO  CALCULATIONS 


The  analytical  methods  described  in  the  previous  chapter  fail  to  account  for 
the  complicated  geometry  involved  in  the  attenuation  of  the  dose  rate  at  positions 
such  as  the  center  and  rear  compartments.  One  method  (and  perhaps  the  only 
comprehensive  one  in  use)  of  including  these  complications  in  the  analysis  is  to 
employ  the  so-called  Monte  Carlo  technique. 

Essentially,  this  technique  involves  the  computation  of  the  probability  of  a 
photon  traversing  a  particular  path  (selected  at  random)  and  striking  the  detector. 
The  results  of  selecting  a  large  number  of  paths  provide  an  estimate  of  the  trans¬ 
mission  probability. 

Such  a  computation  has  been  programmed  by  Raso  and  Woolf®  for  machine 
calculations.  The  geometry  considered  is  an  infinite-plane  slab  of  material  that 
has  photons  incident  on  one  surface,  all  with  given  angles  of  incidence.  The  de¬ 
tector,  located  a  distance  on  the  adjacent  surface  of  the  slab,  counts  only  photons 
entering  it  through  a  solid  angle  given  by  w.  As  seen  in  Figure  16,  this  corre¬ 
sponds  closely  to  counting  only  photons  that  originate  in  a  finite  plane  source  sub¬ 
tending  the  solid  angle  fraction  cj  . 


SOURCE  BEAM 


ANISOTROPIC 

DETECTOR 


Figure  16.  Monte  Carlo  Program  Simulation  of  Slab  Shield  Effects 
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The  results  of  this  first  computEtion  are  directly  applicable  to  measurements 
in  the  front  compartment  of  our  present  egg  crate  experiment.  To  obtain  center- 
compartment  results,  the  program  determines  a  source  of  an  incident-angle  dis¬ 
tribution,  determined  by  the  distribution  leaving  the  first  shield  slab,  and  applies 
this  to  the  second  slab.  In  a  similar  manner,  rear  compartment  results  are  also 
possible .  (A  more  detailed  explanation  of  the  techniques  used  can  be  found  in 
Ref.  8.) 

Results  obtained  by  this  program  are  presented  as 


M<  N 


where  M  is  the  number  of  photons  penetrating  to  the  detector,  N  the  number  of 
photons  incident  on  the  slab  ( "histories  ••),  E,  the  energy  of  the  i  photon  (in  keV), 

A 

and  Qx/p)^  the  mass  absorption  coefficient  of  air  for  the  i  photon.  Now,  if  we 
divide  this  expression  by  (m/p)^E^,  where  the  subscript  refers  to  the  source 
energy,  the  results  are  transformed  to  the  ratio  of  detected  energy  to  incident  en¬ 
ergy  .  This  is  essentially  the  dose  ratio,  and  comparison  may  be  made  with  ex¬ 
perimental  results. 

It  should  be  mentioned  here  that  the  actual  program  used  includes  the  photo¬ 
electric  effect  and  pair  production  as  absorption  processes  only,  thereby  n^lecting 
the  dose  contribution  due  to  pair  annihilation,  bremsstrahlung,  and  X-ray  fluores¬ 
cence.  However,  these  processes  contribute  very  little  to  the  dose  rate,  inasmuch 
as  the  pair  and  photoeffect  cross  sections  represent  only  2%  of  the  total  iron  cross 
section  at  1.25  MeV ,  and  inasmuch  as  these  softer  radiations  are  more  easily 
attenuated  by  the  medium. 

RESULTS  OF  NORMAL  INCIDENCE  COMPUTATIONS 

The  results  of  this  calculation  for  the  centerline  of  the  front  compartment  are 
plotted  in  Figure  17 ,  along  with  the  three  centerline  measurements  for  each  shield 
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FRACTION  OF  BEAM  TRANSMITTED  AS  SCATTERED  RADIATION 


Figure  17 ,  Comparison  of  Monte  Carlo  Prediction  of  Scattered  Dose  Rates 
with  Centerline  Measurements  (Front  compartment) 
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thickness.  Normally  incident  radiation  is  considered.  In  analyzing  this  figure,  it 
should  be  noticed  that  the  origin  is  not  included  in  the  graph,  so  that  differences 
between  Monte  Carlo  and  experimental  results  are  magnified. 

The  20  psf  structure  case  was  the  subject  of  a  further  Monte  Carlo  calculation. 
The  current  angular  distribution,  as  well  as  the  intensity,  was  calculated  at  the 
face  of  the  second  plate,  and  this  result  was  used  as  a  source  to  calculate  dose 
rate  distributions  in  the  center  compartment.  In  addition,  the  current  at  the  face 
of  the  third  plate  was  calculated,  providing  source  for  calculation  of  the  radiation 
in  the  rear  compartment.  At  centerline  detectors,  only  those  dose  rates  due  to 
scattered  radiation  were  obtained.  These  were  presented  in  keV  absorbed  per 
gram  of  detector,  per  second,  and  are  therefore  comparable  to  values  presented 
in  Table  2  through  the  formula 

Dg  =  33.25  -e“^y  (1) 

where  D  is  the  scattered  photon  dose  rate  in  keV/g/sec,  D/D  is  a  dose  ratio 
tabulated  in  Table  2 ,  and  e  ^  is  the  uncollided-photon  portion  of  the  dose  rate . 

3 

The  factor  33.25  equals  10  and  is  introduced  as  described  above.  By 

subtracting  the  uncollided  beam  portion,  a  more  sensitive  comparison  of  the 
structure's  geometry-dependent  portion  of  the  dose  rate  is  possible. 

The  Monte  Carlo  results  are  compared  with  experimental  results  in  Table  16 . 
The  results  for  the  center  and  rear  compartments  neglevt  the  presence  of  side 
walls  and  walls  to  the  rear  of  the  detector.  The  magnitude  of  these  contributions 
may  be  seen  from  the  double  entries  for  the  front  compartment.  Here,  walls  to 
the  rear  of  the  detector  have  been  neglected  for  the  first  entry;  they  are  included 
together  with  side  walls  for  the  second  entry.  The  experimental  values  (with  no 
side  walls)  are  still  higher  than  either  of  the  Monte  Carlo  values.  In  fact,  this 
comparison  is  the  poorest  of  the  three  compartments,  and  was  the  same  compari¬ 
son  that  was  so  poor  in  Figure  17 . 
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TABLE  16 


COMPARISON  OF  MONTE  CARLO  AND  EXPERIMENTALLY  DERIVED 
DOSE  RATES  FOR  NORMALLY  INCIDENT  BEAM  ON  20  PSF  WALLS 


(See  Table  2) 


Compartment 

Detector 

Plane 

(in.) 

! 

Experimental 

Monte  Carlo 

Monte  Carlo, 
plus 

Walls  and  Floors 

2 

8.45 

6.88 

7.22 

Front 

6 

6.78 

6.09 

6.69 

10 

5.12 

4.84 

5.87 

2 

8.21 

7.80 

— 

Center 

6 

6.90 

6.87 

— 

10 

5.22 

5.69 

— 

2 

6.95 

6.81 

— 

Rear 

6 

5.61 

6.13 

— 

10 

4.62 

5.19 

— 

The  value  of  the  uncollided  portion  of  the  dose  rate  in  the  front  compartment 
(with  20  psf  walls)  is  0.58.  Thus ,  an  mcperimental  measurement  of  0 . 84  in  the 
center  of  the  compartment  is  fully  70%  due  to  uncollided  photons .  A  slight  error 
in  the  measurement  results  in  a  large  error  in  the  calculation  of  the  scattered 
portion  in  such  a  situation.  We  believe  this  to  be  the  case  in  the  thin- walled  front 
compartment,  and  thus  prefer  the  Monte  Carlo  result  for  the  reason  explained 
below. 

Consider  a  buildup  factor  of  the  form  (1  +  \)  for  a  beam  of  radiation,  where 
X  is  the  mean  free  paths  of  shield  material.  That  is. 


D  =  D  (1  +  X)  e 
o 


(2) 
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where  D/D^  is  the  ratio  of  the  dose  rate  with  shield  in  place  to  the  dose  rate  with¬ 
out  shield  in  place.  This  form  of  buildup  factor  is  one  of  those  commonly  ac¬ 
cepted,  Since  the  imcollided  photon  contribution  to  the  dose  rate  is  given  by 
e  ^ ,  the  scattered  photon  contribution  is  given  by 

“s  =  •  (3) 

The  value  of  X  that  provides  the  maximum  value  of  D  occurs  when  the  differential 

s 

of  this  expression  with  respect  to  X  vanishes .  Thus , 

(4a) 
(4b) 

In  other  words ,  the  scattered-photon  dose  rate  reaches  its  maximum  1  mfp  into  the 
material,  provided  that  the  buildup  factor  is  reasonably  approximated  by  (1  +  X). 
Since  the  Monte  Carlo  results  in  Table  16  show  this  trend,  they  would  seem  more 
reliable  than  the  experimental  measurements. 

In  the  other  compartments,  the  results  are  more  compatible.  Agreement  is 
within  10%,  and  in  most  cases  within  5%.  Here,  it  is  difficult  to  assign  these  small 
deviations  to  underlying  difficulties  with  either  method,  and  perhaps  it  would  be  a 
useless  task  to  attempt  to  do  so. 

RESULTS  OF  SLANT  INCIDENCE  COMPUTATIONS 

The  results  of  the  angular  runs  (Tables  8  through  12)  were  compared  with 
Monte  Carlo  calculations  for  the  centerline  detectors  of  the  front  compartment. 
These  comparisons  are  presented  in  Table  17 .  The  Monte  Carlo  code  presumed  a 
constant  intensity ,  but  since  we  reduced  the  intensity  of  radiation  incident  on  the 
front  plate  by  tilting  the  experimental  structure,  the  appropriate  form  of  Eq.  (1) 
for  these  comparisons  of  a  40  psf  structure  is 


e"^(l  -  X)  =  0 
X  =  1  . 
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TABLE  17 


COMPARISON  OF  MONTE  CARLO  AND  EXPERIMENTALLY  DERIVED 
DOSE  RATES  FOR  SLANT  RADIATION  ON  40  PSF  STRUCTURE 


COS  0 

Detector  Plane  (in.) 

2 

6 

10 

Monte  Carlo 

1.0 

7.91 

6.94 

5.72 

0.75 

8.38 

6.68 

4.71 

0.50 

8.47 

6.03 

3.66 

0.25 

6.25 

4.15 

2.36 

0.00 

2.08 

1.59 

0.99 

Experimental 


1.0 

8.99 

7.33. 

6.00 

0.96593 

8.72 

7.34 

6.65 

0.86603 

9.32 

7.40 

7.01 

0.7071 

9.55 

7.66 

2.02 

0.5000 

-  —  -  —  — 1 

9.62 

5.97 

2.52 

In  general,  it  may  be  seen  that  agreement  is  quite  good.  The  40  psf  walls 
correspond  to  approximately  1  mfp,  which  reduces  the  experimental  difficulty 
discussed  earlier  in  this  chapter  under  normal  incidence  computations ,  and  does 
not  provide  too  great  a  difficulty  for  the  code  in  calculating  histories . 

ISOTROPICALLY  INCIDENT  RADIATION 

Combining  the  results  of  several  Monte  Carlo  calculations  for  centerline  dose 
rates  due  to  radiation  incident  at  various  angles  gives  a  reasonable  estimate  of  the 
total  dose  rate  due  to  isotropic  radiation.  (A  similar  effort  was  described  in 
Chapter  4,  for  the  experimental  results.)  Only  the  scattered  photon  component  was 
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considered  in  the  computations,  since  the  direct  component  may  be  easily  com¬ 
puted  analytically,  as  was  done  above. 

So  far  in  this  chapter  we  have  considered  only  parallel  solutions  for  a  partic¬ 
ular  problem  by  using  both  Monte  Carlo  computations  and  experimental  measure¬ 
ments  .  Monte  Carlo  computations  were  performed  for  isotropically  incident 
radiation  for  structures  of  20  psf  and  40  psf  walls .  It  is  felt  that  the  previous  sec¬ 
tions  justify  a  measure  of  reliance  in  the  agreement  of  results  of  the  two  methods. 
Therefore,  the  two  calculations  for  structures  of  differing  wall  thicknesses  should 
present  a  limited  range  of  evaluations  of  the  isotropic  radiation  problem.  How¬ 
ever,  because  of  the  ease  of  solution,  the  isotropic  radiation  response  was  com¬ 
puted  for  40  psf  walls  for  the  front  compartment. 

Figures  18  through  21  present  the  constructed  isotropic-source  reduction  factor 
versus  the  solid-angle  fraction  co  for  the  various  wall  thicknesses  and  compartments 
of  the  structure .  The  Monte  Carlo  and  experimental  points  are  compared  in  each 
figure  with  the  corresponding  transport  solution  of  Ref.  7.  This  affords  a  two-way 
comparison  of  data  with  the  basic  curves  for  the  roof  source  reduction  factor  of  the 
OCD  Engineering  Manual.  To  compare  these  data  with  Spencer's  transport  solu¬ 
tion,  it  is  necessary  to  subtract  from  his  solution  LpC)L  (X,  w)— the  component  due 
to  uncollided  radiation.  This  leads  to  the  equation 


Dg  =  BE^(X^)L(X)L^(X,  CO)  - 


E^(X)  -  E^(X/l-(0) 


(6) 


Here,  X  is  the  shield  thickness  in  mean  free  paths,  and  is  therefore  proportional  to 

X.  The  exponential  function  Ej^(X)  enters  the  expression  as  it  did  in  the  previous 

section.  The  first  term  on  the  right  in  Fq.  (6)  represents  the  total  dose  rate,  and 

the  second  represents  the  unscattered  photon  contribution.  The  dose  rate  in  air  3ft 

from  an  infinite  plane  of  this  source  is  given  by  the  expression  BEj^ ,  (X^) ,  where 

X^  is  the  shield  thickness  of  the  air  under  standard  conditions.  The  function  D 
o  s 

defined  by  Eq.  (6)  is  plotted  as  the  transport  solution  in  Figures  18  through  21.  We 

12 

obtain  the  buildup  factor  B  from  a  note  by  Berger  and  Spencer.  Thus,  except  for 

the  exponential  integrals,  the  expression  for  D  is  derived  from  transport  calculations. 
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CHAPTER  6 


CONCLUSIONS  AND  RECOMMENDATIONS 


CONCLUSIONS 


SHIELDING 

Several  general  conclusions  may  *'e  drawn  regarding  the  gamma  shielding 
offered  by  a  compartmented  structure. 

1 .  It  appears  evident  that  backscatter  represents  a  small  contri¬ 
bution  to  the  dose  rate.  This  was  on  the  order  of  2%  in  the 
experiment ,  although  in  a  situation  where  a  structure  composed 
of  a  lower  atomic  number  material  is  exposed  to  somewhat 
lower  energy  quanta  of  radiation,  the  contribution  would  be  ex¬ 
pected  to  be  somewhat  higher. 

2.  The  experimental  results  indicate  that  dose  rates  are  dependent 
not  only  on  the  shield  thickness  in  front  of  a  given  detector ,  but 
also  on  the  distribution  of  the  shield.  Detectors  in  each  com- 
partme .'t,  behind  the  same  mass  thickness,  show  a  marked 
drop  in  dose  rate  as  the  distance  from  the  shield  slab  increases . 
Further,  the  experimental  data  {Tables  2  through  4)  show  one 
shield  slab  is  generally  less  effective  than  two  slabs ,  each  half 
as  thick. 

METHODS  OF  SOLUTION 

The  figures  and  tables  in  thie  report  show  that  there  is  good  agreement  among 
the  three  methods  of  solution- experimental,  Monte  Carlo,  and  moments  solution 
of  transport  theory.  Thus,  any  problem  of  gamma  shielditig  by  thin  shields  should 
be  solvable  by  any  of  the  three  methods,  given  sufficient  time  and  care,  for  there 
are  no  surprises  in  the  physics  involved.  However,  differences  arise  when  one 
considers  the  effort  required. 

1 .  The  transport  solution  probably  takes  the  most  effort ,  but  gives 
answers  that  are  more  self-consistent,  if  not  truer.  Difficulties 
are  encountered  only  when  the  context  ui  the  problem  is  changed; 
for  example,  we  did  not  have  such  a  solution  to  compare  with  ttie 


46  •  u 


R  L  t  ■  ■  T  O  N 


HAtlArHUflT 


effect  of  looking  at  a  finite  cross  section  of  a  shielded  mono- 
directional  source.  Such  a  case  would  require  an  extension 
of  the  graph  collection  nf  the  moments  method  solutions. 

2.  The  Monte  Carlo  program  appears  to  be  excellent  for  its  ir- 
tended  purpose  — the  solution  of  a  compartmented  structure 
shielding  effect.  Again,  however,  the  method  requires  time 
and  patience  and,  therefore,  far  from  all  of  the  experimentally 
measured  points  were  checked  by  the  computer  calculations. 

3.  In  the  ultimate  analysis,  the  experimental  approach  is  the  best, 
even  though  the  experiment  here  is  far  from  ideal.  It  is  diffi¬ 
cult  to  assess  the  effects  of  humidity,  background,  and  mech- 
ical  shock  on  each  dosimeter's  residual  charge.  The  Styrofoam 
holders  present  a  small,  variable  shielding  effeet.  Finally, 
small  systematic  errors  in  the  measurement  of  a  total  exposure 
can  lead  to  gross  errors  in  the  calculated  exposure  to  scattered 
photons ,  particularly  when  the  unscattered  component  is  much 
greater  than  the  scattered.  As  a  result,  the  internal  consist¬ 
ency  of  the  experimental  results  is  considerably  lower  than  in 
the  other  cases.  Nevertheless,  the  versatility  of  the  experi¬ 
mental  approach,  its  independence  of  the  mathematical  model 
that  must  be  otherwise  developed,  and  Its  economical  aspects 
make  it  a  formidable  tool  for  cross-checking  a  previous  solu¬ 
tion  as  well  as  researching  an  unknown  one . 

RECOMMENDATIONS 

1.  The  Monte  Carlo  code  for  attenuation  of  radiation  in  a  compart¬ 
mented  structure  appears  to  be  of  an  accuracy  to  warrant  its 
further  use  in  similar  shielding  studies  of  compartmented 
structures . 
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2.  This  study  indicates  that  the  OCD  Engineering  Manual  over¬ 
estimates  the  dose  rate  obtained  behind  multiple  thin  layers. 

An  experiment  with  isotropic  sources  on  the  face  of  the  shield, 
with  a  variable  number  of  interposing  layers,  should  be  per¬ 
formed  to  further  explore  this  discrepancy.  An  analysis  should 
be  made  to  see  if  the  rooftop-source  reduction  factor  C^,  de¬ 
fined  in  the  Engineering  Manual,  need  be  modified  to  account 

for  interior  floor  spacing  in  a  multistoried  building.  (Block- 
13 

house  experiments  have  indicated  that  the  theoretical  reduc¬ 
tion  factor  provides  quite  a  conservative  estimate  of  the  experi¬ 
mental  factor  when  floors  exist  between  roof  and  detector . ) 
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